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Abstract 

We consider a two-Higgs doublet scenario containing three SU(2)l singlet heavy neutrinos 
with Major ana masses. The second scalar doublet as well as the neutrinos are odd under a Z 2 
symmetry. This scenario not only generates Majorana masses for the light neutrinos radiatively but 
also makes the lighter of the neutral Z 2 -odd scalars an eligible dark matter candidate, in addition to 
triggering leptogenesis at the scale of the heavy neutrino masses. Taking two representative values 
of this mass scale, we identify the allowed regions of the parameter space of the model, which are 
consistent with all dark matter constraints. At the same time, the running of quartic couplings 
in the scalar potential to high scales is studied, thus subjecting the regions consistent with dark 
matter constraints to further requirements of vacuum stability, perturbativity and unitarity. It is 
found that part of the parameter space is consistent with all of these requirements all the way up to 
the Planck scale, and also yields the correct signal strength in the diphoton channel for the scalar 
observed at the Large Hadron Collider. 


1 Introduction 


The discovery of c a Higgs-like boson’ [1,2], apparently completes the Standard electroweak model 
(SM). However, it also compels us to mull over the necessity of physics beyond the SM. To mention 
one motivation for this, a SM Higgs of mass around 125 GeV, can break the absolute stability of vacuum 
at 10 8-9 GeV, if the top quark mass ( M t ) and the strong coupling constant (<a s ) are on the upper sides 
of their respective uncertainty bands. A recent next-to-next-to-leading order (NNLO) study [3] reveals 
that absolute stability up to the Planck scale requires 


M h [GeV] > 129.4 + 1.4( M ‘l Qe yl- 17;U ) - ± 


0.7 


0.0007 


( 1 ) 


which, with all theoretical and experimental errors in Mt and a s , implies that the absolute vacuum 
stability of SM is excluded at 95% CL for M\ < 126 GeV. To compensate the negative contribution 
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from the fermion loops (mainly the top quark loop), which rapidly drives the quartic self-coupling of 
Higgs to below zero, we need extra scalar loops that can ameliorate the problem. This provides one of 
the important motivations to consider a theory with new scalars which can make the vacuum stable 
up to high scales while having consistent results at low scales. 

At the same time it is worthwhile to remember two rather pressing issues which prompt one to look 
beyond the Standard Model (BSM). These are the non-zero mass and mixing of neutrinos and the 
likely existence of Weakly Interacting Massive Particle(WIMP), contributing to Dark Matter (DM). 
While various Solar, atmospheric, reactor and accelerator-based experiments provide strong evidence 
for non-zero masses of neutrinos and mixing among their different flavors, several astrophysical and 
cosmological observations suggest the existence of some exotic particles that constitute the DM content 
of the universe. Given all this, it is an attractive idea to look for a theory that can simultaneously 
address all of the aforesaid problems in one framework. Here we consider one such scenario. 

We investigate a model, first proposed in reference [4], that extends the SM with an extra Higgs doublet 
and three right handed neutrinos with a Z 2 symmetry, under which all SM particles are even while 
this additional scalar doublet and the right-handed neutrinos are odd. This symmetry prevents the 
additional doublet from having a vacuum expectation value(vev) thus vetoing the tree-level neutrino 
mass generation. Moreover, there exists a stable scalar particle in the form of the lightest neutral 
mass eigenstate of the additional doublet, which yields an eligible DM candidate. The extra doublet is 
essentially an Inert Doublet. Although a lot of study has already taken place on minimal Inert Doublet 
models [5-23], the extra appeal of this model lies in the radiative generation of neutrino mass. Though 
various aspects of this scenario have already been investigated [24-35], the present study includes the 
following points which have not been emphasized before. 

• We investigate the vacuum stability of this model at various scales, and identify the regions of its 
parameter space, which keeps the model valid all the way unto the Planck scale. The contribution 
of additional scalar fields as well as the right-handed neutrinos to the renormalisation group (RG) 
equations [given in Appendix A] has been taken into account here. Using these modified RG 
equations we evaluate the scalar quartic couplings at different scales. During the evolution of the 
quartic couplings we demand not only vacuum stability but also perturbativity of the couplings 
as well as unitarity of the 2 —)> 2 scattering matrix at each scale. 

• The heavy right-handed neutrinos introduce a new mass scale (M) to the theory and the neutrino 
Yukawa couplings contribute to the RG evolution only at this mass scale. This brings out greater 
implications on the parameter space that distinguishes this model from minimal Inert doublet 
models. In this paper, we will show the salient features of the model that emerges from the above 
fact. 

• This model also contributes to the study of leptogenesis due to this heavy right-handed fermions. 
The values of the right-handed neutrino mass scale used in the high-scale analysis are taken to be 
commensurate with leptogenesis constraints. Thus the part of the model space consistent unto 
the Planck scale is also supportive of leptogenesis. 

• We also examine the candidature of the lightest Z 2 -odd particle as DM candidate, and identify 
the allowed values of the couplings yield the right relic abundance. We ensure that the DM 
candidate is consistent with the recent result of direct detection experiments. We identify a 
substantial region of the parameter space, which simultaneously satisfy the vacuum stability, 
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perturbativity and unitarity requirements and accommodate the DM candidate with the correct 
relic density. 

• Finally, we examine the 125 GeV scalar and make sure that the signal strengths in the observed 
channels (such as diphotons) are consistent with data from the Large Hadron Collider (LHC). 

We organize the paper as follows. In Section 2, we briefly describe the model and its various features. 
In Section 3, we explain all the theoretical constraints and collider constraints that we use in the RG 
running of different quartic couplings. Next, in Section 4, we discuss the DM aspects of this model. 
After explaining our analysis strategy in Section 5, we present our results related to high-scale validity 
in Section 6. Finally in Section 7, we summarise our results. 


2 The Radiative Neutrino Mass Model with an Inert Doublet 


In addition to the SM fields, the radiative neutrino mass model with an inert doublet [4], contains 
a Higgs doublet ($ 2 ) and three right handed (RH) neutrinos (A^) with an unbroken Z 2 symmetry, 
under which the doublet and the right handed neutrinos are odd while all other SM particles are even. 
Being odd under the symmetry, $2 does not acquire any vacuum expectation value (vev) and has no 
tree-level couplings to fermions. 

The relevant Yukawa and mass terms are 

-Cy = y ij mie j + h.c+^(N?N i + h.c),{i,j = 1-3) (2) 


Here £i are the left-handed lepton doublets and Mi are the Majorana mass term for the heavy right- 
handed neutrinos Ni 
The scalar potential is 


V = 
+ 


y(<^2) 2 


+ A2 ($2^2 
+ h.c. 


+ A 3 ($I$i)(^$2) + A 4 (^$i)($I$ 2 ) 




( 3 ) 


where all parameters are real, and 4>i is the SM Higgs doublet. 

The two scalar doublets can be written as 

$i= ( “ d ' *’■-(*(£*>) 


( 4 ) 


where, v = 246 GeV, is the electroweak vacuum expectation value (vev). One thus has five physical 
states (A, iL, A, H^) and three Goldstone bosons (G, G±). While h corresponds to the physical SM- 
like Higgs field, the inert doublet components are one CP-even neutral scalar(iL), one CP-odd neutral 
scalar(A) and a pair of charged scalars (H±). The physical masses are given by 

M 2 h± = m| 2 + L 3 v 2 

M‘h = m % 2 + 


3 




( 5 ) 


M\ = m% 2 + A av 2 

where A l /a = \ (A 3 + A 4 ± A 5 ). The value of Ai is determined using M/, = 125 GeV. 

Majorana masses for the light neutrinos are generated radiatively through one-loop exchange of the 
Z 2 -odd neutral scalars. The general expression for the loop-induced contributions to the light neutrino 
mass matrix [4] is 

kav A UikVjkMk 

ij ^ 16vr 2 

k= 1 

Thus, the neutrino masses and mixing are determined by the inert scalar masses and the right-handed 
neutrino masses Mi. These masses represent the scale of lepton number violation and hence that of 
leptogenesis [36,37] in this model. Some studies have already been done in this context [11,30,32,33, 
38]. Our choice of the right-handed neutrino (Majorana) mass scales made our study automatically 
compatible with leptogenesis. To satisfy the necessary constraints in the low DM mass region where 
Mdm < Mw, one must take the lightest of the Majorana masses to be M > 110 TeV, where as in 
the high DM mass region (M^m > 500 GeV) the bound is only 1 TeV [11]. Hence, to be consistent 
in both cases, we use two values of M, (a) M — 110 TeV and (b) M — 10 9 TeV. While, in one hand, 
choice of (a) is motivated by the idea of having the lowest possible leptogenesis scale, we choose to 
work with (b) which have interesting consequences on the RG runnning. We will show in later sections, 
how the mass scale (b) of M affects the stability of the vacuum mainly in the high DM mass region 
and eventually explain the physical reasons behind it. 

Along with the above restrictions we also demand M v ^ 0(0.1 eV) to be consistent with neutrino 
oscillation data for some fixed M and other exotic scalar masses. However, for simplicity, we consider 
only one diagonal Yukawa coupling (y„) and do not look into the hierarchical details of the Yukawa 
matrix. At this point, it is to be noted that a more rigorous study with the intricate flavor structure of 
the neutrino Yukawa matrix can highlight the region of the model space that fits the observed pattern of 
neutrino mixing [35]. However, we would like to emphasize that the broad conclusions on the high-scale 
validity of this scenario vis-a-vis the DM constraints remain unchanged. Finally, we should mention 
that the lighter state between H and A is the DM candidate. We present our illustrative results for 
cases where H plays this role. 

3 Constraints from perturbativity, unitarity, vacuum stability and 
collider data 


Mfi 


In 


Mi 


H 


Ml 


In 


Mi 


Ml-Ml Ml M 2 A -Ml Ml, 


( 6 ) 


In this section, we briefly describe the constraints that are imposed on the model parameters and how 
exactly they shape the results so obtained. 


3.1 Vacuum stability 

The scalar potential is considered bounded from below, if it does not turn negative for large field values 
along all possible field directions. In this model, stability of the electroweak vacuum is ensured up to 
some specified energy scale if the following conditions are satisfied for all scales Q up to that scale: 

vscl : Ai(Q) > 0 (7a) 
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vsc2 : 

HQ) > 0 

(7b) 

vsc3 : 

A3 (Q) + >/M(Q)^ 2 (Q) > 0 

(7c) 

vsc4 : 

Aa(Q) + A 4 (Q) — |A 5 (Q)| + VAi(Q)A 2 (Q) > 0 

(7d) 


Such conditions have been elaborately discussed in literature [39-42]. One should make a note that 
these conditions ensure absolute stability of the electroweak vacuum. For metastability, the conditions 
are somewhat less stringent. 


3.2 Pert ur bat i vity 

For the scalar quartic coupling A^(i = 1 — 5), the criterion for perturbativity is, 

A i(Q) < 47r (8) 

The corresponding constraints for the Yukawa and gauge interactions are, 

Vi(Q), 9 i(Q ) < (9) 

where, Q represents the energy scale at which they are being computed. We demand that the criteria 
in Eq. 8 be obeyed at all energy scales up to the cut-off of this model. 

3.3 Unitarity 

A further set of conditions come on demanding unitarity of the scattering matrix comprising all 2 —» 2 
channels involving, by the optical theorem [43-46]. In our context, this translates into the condition 
that each distinct eigenvalue of the aforementioned amplitude matrix be bounded above at 8tt (after 
factoring out from the matrix). These eigenvalues are: 


a± = -(Ai + A2) ± \j -(Ai — A2) 2 + (2A3 + A4) 2 
b ± = i(A 1 + A 2 )±^(A 1 -A 2 )2 + Al, 

c± = d ± = \( Ai + A 2 )±yt (Al _ A2) 2 + A 2 ) 

e i = (A3 + 2X4 — 3A5), 
e2 = (A3 — A5), 

fl = /2 = (V + A4), 

/+ = (A3 + 2A4 + 3A5), 

/- = (A 3 + A5). ( 10 ) 


3.4 Collider data 

In addition to the theoretical constraints discussed above, important bounds on scalar mass parameters 
come from collider data. 
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• In order to identify h with the SM-like Higgs as observed by the ATLAS and CMS collaborations, 
one requires M\ ~ 125 GeV. 


• To be consistent with the LEP bounds, one must have 


Mh + Ma > Mz , ( 11 ) 

M h ± + Ma/h > Mw • 

Moreover, for neutralino as in the supersymmetric context, LEP-II searches limit pseudo-scalar 
mass (Ma) to 100 GeV when Mh < Ma [8,47]. Similarly, chargino search data, properly 
extrapolated, imply M h ± > 70 GeV [48]. 


• Though all the tree-level couplings of h are identical to those of the SM Higgs, the charged scalar 
H ± potentially modifies the loop induced couplings /177 and / 1 Z 7 via loop contributions [49-53]. 
We theoretically compute the signal strength /i 77 for h decaying to the diphoton channel as the 
ratio of the decay width in the ‘model’ to that in the SM. Its experimental value reported by 
the ATLAS and CMS collaborations stand at 1.17 ± 0.27 and 1.13 ± 0.24 respectively [54,55]. 
Demanding the signal strength to be within 2 a limits of the experimentally quoted values puts 
further constraints on the model. We use the limit on /i 77 weighted as below: 


-7 = (— 9 ) ATLAS + (-n)cMS 
cr z 



( 12 ) 

(13) 


where, the numerators on the right hand side denotes the central values of the respective exper¬ 
imental results and a are the corresponding uncertainties. 


4 Dark Matter Issues 


As stated earlier, we identify H as the DM candidate. For complimentary choice, namely, A with the 
same mass as the DM candidate, we have checked that the contribution to the relic density is of similar 
magnitude. The relevant DM constraints to be considered are as follows: 

• According to recent PLANCK experiment [56] the observed cold DM relic density is 

n DM h 2 = 0.1199 ± 0.0027 (14) 

We restrict our result up to 3cr deviation from the central value. 

• Strong constraints come from direct DM searches. Recently, XENON100 [57] and LUX [58] 
experiments have put upper-bound on the DM-nucleon scattering cross-section for a wide range 
of the DM mass. In our case, the direct detection cross-section strategy is based on t-channel 
Higgs mediation. We choose to work in the region of the parameter space allowed by the LUX 
limit. 

• For Mh < Mh/ 2, the decay mode of Higgs to two DM particle (h —» HH) will presumably 
contribute to the invisible decay width of the Higgs boson. We take into account the current 
constraint on the Higgs invisible branching ratio (< 20%) from model independent Higgs precision 
analysis [59,60]. 
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• Our parameter scan has included situations where the Z 2 ~odd scalar and the pseudo-scalar are 
close by in mass and with co-annihilation yield correct relic abundance within the ±3 a of the 
observed values when the DM mass is high (> 500 GeV). For DM mass in between 100 GeV 
and 500 GeV, there is prima facie the possibility of co-annihilation between the spinless DM 
candidate and a right-handed neutrino which can give rise to correct relic abundance, as discussed 
in Ref. [61]. In our case, however, such a possibility is ruled out by the additional requirement of 
leptogenesis, which in turn puts a lower bound on the mass range for the right-handed neutrino(s). 

Hence, only the following two cases survive all the above constraints. 


4.1 Case-A: Low mass DM (50 GeV < M H < 90 GeV) 

In this mass region, the dominant annihilation channel for H goes to the DM self-annihilation processes 
through h mediation. This keeps the relic density at the right level. For both positive and negative 
values of DM-DM-Higgs couplings, the relic density remains within the allowed range as long as the 
Mh < 90 GeV. The sub-dominant contribution to the relic density comes from the t-channel processes 
to vector boson final states mediated by A and H ± . A detailed discussion in this regard on a similar 
model can be found, for example in [11,13,22]. However, the coupling A 2 has no effect in the relic 
density calculation. In the next section, we will discuss this results elaborately. One more notable 
point is that, annihilation processes that mediated by the heavy right-handed neutrinos give negligible 
contribution (less than 1%) to the relic density calculation. These processes are suppressed by the 
heavy mediator mass. 


4.2 Case-B: High mass DM ( M H > 500 GeV) 

The interesting feature of this region is that, the correct relic abundance can be achieved if and only 
if H, A and are almost degenerate, at most have a mass difference of the order of 10 GeV. This 
is mainly because at this high mass, the annihilation channels with vector boson final states can have 
direct quartic couplings (. HHVV , V = W ± , Z ) or can be mediated by any of the three scalars through 
t/u channels. These diagrams yield too large an annihilation cross-section to match with the proper 
relic density. However, cancellation between direct quartic coupling diagrams and t/u channels diagram 
occurs for mass-degenerate i7, A and H ± 1 which in turn brings down the annihilation cross-section to 
the desired range. 


5 Analysis strategy 


The aim of this study is to probe the parameter space of an inert doublet model (IDM) augmented 
with heavy RH neutrinos compatible with various theoretical and experimental constraints elaborated 
in the previous sections. We carry out our investigation in the two separate mass regions. In each 
region, we scan over the relevant parameters, namely, the masses M#, Ma and M^, and the coupling 
A l. With Mh fixed at 125 GeV, A 2 , A£, Mh , Ma and M^ fix all the remaining quartic interactions. 
At this point it is to be mentioned that the parameter A 2 can not be constrained by any physical 
observable at least at tree level. However, A 2 is bounded by the stability and perturbativity condition 
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and we have explicitly checked that 0 < A 2 < 0.36 to satisfy the theoretical constraints. Rather than 
scanning over A 2 , we have demonstrated our results with two benchmark values (0.1 and 0.001) for it, 
both of which are well within the above limit. 

These quartic couplings are then used as the electroweak boundary conditions at Q — M t and their RG 
evolution to high scales is studied. Here M t denotes the top quark pole mass. The reader is reminded 
that the effect of the RH neutrinos is turned on at a scale Q = M. Thus, for M t < Q < M, we do not 
include the RH neutrino contributions to the one-loop RG equations. We include such conditions for 
Q > M and use the RG equations listed in the appendix. The scale up to which the scenario remains 
consistent is denoted by A uv- For a generic parameter point A i(Q — M^), we check the aforementioned 
theoretical constraints at all intermediate scales up to A jjy. If the constraints are all satisfied, we tag 
A i(Q — M t ) as an allowed point. This marks out an allowed region in the parameter space defined 
at the electroweak scale. Moreover, the effects of constraints stemming from the DM observables and 
collider searches are examined independently in this region. The finally allowed parameter regions 
are thus identified and presented for benchmark values of A 2 and M. We use the publicly available 
package micrOMEGAs-v3.6.9.2 [62] for DM analysis. In all the cases, we also make sure that our 
allowed parameter space does not violate the oblique T-parameter constraints [63,64]. 

Amongst the SM fermions, only the top quark plays the dominant role in the evolution of the couplings. 
The boundary condition for its Yukawa interaction at the electroweak scale is fixed by y t (M t ) = 
^ Mt (1 — a s (M t )). We have used M t = 173.39 GeV throughout our analysis. 

6 Numerical results 

6.1 Low mass DM region 


We perform a detailed parameter space scan where Mh < 100 GeV. In this scan, 
bounds as discussed in Sec 3.4. 

we impose the LEP 

X L : [-1.0,1.0] 

(15) 

M h : [50.0 GeV, 90. GeV] 

(16) 

M a : [100. GeV, 500. GeV] 

(17) 

M h + : [100. GeV, 500. GeV] 

(18) 


We solve the RG equations for two values of M, 110 TeV and 10 9 TeV respectively. We then show the 
allowed parameter space in the A l~ Mh plane for different choices of A jjv in Fig 1 and 2. The regions 
denoted by A (red), B (cyan) and C (green) denote A jjy = 10 6 , 10 16 and 10 19 GeV respectively. We 
overlay the region allowed by the Higgs to diphoton signal strength within 2a limits of the current 
data as grey region named D. As mentioned earlier, the full analysis is done for two values of A 2 at the 
electroweak scale (0.1 and 0.001). 

Let us briefly summarize the features that emerges from the figures and their possible explanations. 

• From Fig.l, it can be seen that an Inert doublet model with right handed neutrinos posses stable 
vacuum even up to the Planck scale, but the higher the cut-off scale A uv, the less amount of 
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Figure 1: Regions compatible with the theoretical constraints in the Mh~\l plane for M — 110 TeV 
(left panel) and 10 9 TeV (right panel) with three different choices of A uv and two values of A 2 (upper 
and lower panel) . The regions denoted by A (red), B (cyan) and C (green) obey these constraints 
up to A uv = 10 6 , 10 16 and 10 19 GeV, respectively. The grey region denoted by D keeps the Higgs to 
diphoton signal strength within 2 a limits of the current data. 


parameter space becomes allowed. To understand this, one can recollect that the top quark 
Yukawa coupling in the SM is responsible for the downward evolution of the scalar self-coupling, 
which poses a threat to the vacuum stability. The presence of the additional scalar quartic 
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couplings (A 's) in a model like this offsets such an effect; however, the boost thus provided to 
these couplings tend to violate the perturbativity and unitarity condition. This necessitates a 
tightrope walking, and the scale up to which it is possible is Ajjv- Hence, it is natural that for 
higher AjjVi fewer combinations of parameters will achieve this fine balance, and consequently 
the allowed region shrinks. 

• Fig.l also shows that A^ is bounded on both sides and the limits stay almost same for different 

right handed neutrino mass scale(M). This is of no surprise and can be easily understood. Since 
A 3 = —Mjj+ A^u 2 ), the upper bound on A^ is imposed by the requirement of perturbative 

unitarity. This is because a higher positive value of A^ makes A 3 large at the electroweak scale 
which violates perturbative unitarity in the course of its evolution under RG. On the other hand, 
a large negative value of A^ induces a large negative value to A 3 . As a consequence, the vacuum 
stability condition vsc4 of Eq.7d is violated even near the electroweak(EW) scale. This puts 
a lower limit on A^ independent of M and Ajjvi as evident from Fig.l. However, it must be 
noted that the lower limit of A^ is not independent of A 2 , which is another consequence of Eq.7d 
that requires the condition A 3 + A 4 — |A 5 1 > — \A 1 A 2 to be satisfied. With the decrease in the 
value of A 2 the lower limit of the combination A 3 + A 4 — | A 5 1 enhances. This combination can be 
approximated to A^ in the perturbative regime and hence lower values of A 2 increases the lower 
bound of Al, as can be seen from the upper and lower panel of Fig.l that corresponds to A 2 = 0.1 
and A 2 = 0.001 respectively. 

• To check the compatibility of the DM constraints with the theoretical ones, we look for the 
region allowed by the 3cr limits on Ddm ^ 2 from PLANCK data and 90% CL limit on the spin- 
independent DM-nucleon scattering cross-section from LUX data. In Fig. 2, we show the param¬ 
eter space allowed by the entire set of DM constraints. An inspection of Fig. 1 and Fig. 2 shows 
that almost the full parameter space allowed by the DM constraints lies in the region which is 
also favoured by the vacuum stability condition all the way up to the Planck scale. However, the 
region corresponding to A^ < —0.1 (shaded region in Fig. 2) does not lead to the stable vacuum. 

• It should be noted that in Fig. 2 we keep M^ and Ma fixed at 200 GeV. For higher values of 
M h ± and Ma , there is hardly any change in the annihilation cross-section. However, for values of 

and Ma less than 200 GeV, the allowed region of Fig. 2 gets slightly modified. For example, 
for Mh — 70 GeV and M H ± = Ma = 200 GeV, the DM-DM-Higgs coupling Xl — 0.007, but for 
M h ± — Ma — 100 GeV, one needs Xl — 0.009 to satisfy the relic density constraint. It should 
be noted however that both of the above points in the parameter space are within the stability 
region as shown in Fig. 1. 

Therefore, it is not possible to constrain M^ and Ma using DM constraints alone, theoretical 
constraints however predict strong upper bounds on these masses, as is evident from Fig. 3 and 
Fig. 4. 

• In Fig. 3 and Fig. 4, we show the theoretically valid regions in the A 3 — M^ and A l~ Ma planes. 
As in previous cases, we exhibit the same for two different M and A 2 values. Also, in each case, 
we overlay the parameter spaces allowed by the 2a limit of Higgs to diphoton signal strength 
(region D). We observe that a tight upper bound of 160-170 GeV is realised on the masses M^ 
and Ma for the cut-off at the Planck scale and the couplings (A 3 , A l) are also bounded. The 
upper bounds on the masses M^ and Ma are imposed by the requirement of perturbativity and 
unitarity up to the desired cut-off. With Mh in the aforementioned range, large masses of the 
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Figure 2 : Region allowed by imposing the constraints on relic density(RC) and spin-independent 
cross-section(SI) for DM-nucleon scattering. The red(gray) region is allowed only by the requirement 
of f^DM ^ 2 being in the correct range. The black region is allowed by both the Ddm ^ 2 and direct 
detection constraints. The shaded horizontal band below is disallowed by vacuum stability conditions. 
Here, M H ± — M A = 200 GeV. 


other Z 2 -odd scalars imply high values of the quartic couplings at the electroweak scale which 
potentially violate perturbativity or unitarity at some high scale. On the other hand, the vacuum 
stability condition A 3 + \A 1 A 2 > 0 forbids large negative values of A 3 and perturbativity puts 
the upper bound. The explanations of the limits of A^ have already been mentioned above. 

• The new physics contribution to the Higgs to diphoton decay channel comes from the charged 
Higgs loop which is a function of M H ± and A 3 . The fact that A 3 can not be large negative 
results in a fall in the signal strength for diphoton channel. However, the parameter space for 
A jjv = 10 19 GeV is still allowed by the current limits on /i 77 defined in Eq. 13. 

6.2 High mass DM region 

This section demonstrates the high scale validity of our scenario in the limit of a high DM mass. As 
discussed earlier, one needs to tune the mass splitting amongst i7, H ± and A and the coupling A^ to 
an appropriate degree in order to achieve a relic density within the desired bounds. It is seen that the 
maximal allowed splitting (AM) amongst the masses of the Z 2 odd scalars is 10 GeV. As previously 
mentioned, for each chosen values of A 2 (0.1 and 0.001), one is thus motivated to scan the high DM 
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Figure 3: Regions allowed by the theoretical constraints projected in the M^- A 3 plane for M — 110 TeV 
(left panel) and 10 9 TeV (right panel) with three different choices of A jjv and two values of A 2 (upper 
and lower panel). The regions denoted by A (red), B (cyan) and C (green) obey those constraints up 
to A uv = 10 6 , 10 16 and 10 19 GeV, respectively. The grey region denoted by D shows the 2 a allowed 
limit of the Higgs to diphoton signal strength. 
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Figure 4: Regions allowed by the theoretical constraints projected in the Ma~Xl plane for M — 110 TeV 
(left panel) and 10 9 TeV (right panel) with three different choices of A jjv and two values of A 2 (upper 
and lower panel). The regions denoted by A (red), B (cyan) and C (green) obey those constraints up 
to A uv = 10 6 , 10 16 and 10 19 GeV, respectively. The grey region denoted by D shows the 2 a allowed 
limit of the Higgs to diphoton signal strength. 
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mass region in the following ranges: 


A l : [- 1 . 0 , 1 . 0 ] 

(19) 

M h : [550.0 GeV, 1000.0 GeV] 

( 20 ) 

Ma • [Mu, Mh + 10.0 GeV] 

( 21 ) 

M]j+ : [Mh, Mh + 10.0 GeV] 

( 22 ) 


Unlike the previous case, while the theoretical constraints ruled out a large portion of the parameter 
space, the DM constraints put a less stringent bound on it in this high DM mass region. Therefore, in 
Fig. 5 and Fig. 6 , we demonstrate the parameter region allowed only by the DM constraints and then 
overlay the part which are further allowed by theoretical constraints. The full (pink) region denoted 
by RC+ SI shows the valid parameter space allowed by DM constraints. In accordance with previous 
notation, the regions denoted by A (red), B (cyan) and C (green) denote A uv = 10 6 , 10 16 and 10 19 
GeV respectively. 

Let us explain the various features of the model that emerges from the figures, in detail. 

( 1 ) The interplay of the theoretical and experimental constraints is studied in the form of correlation 
plots in the Mh -A l (Fig. 5) as well as M H ± - A 3 (Fig. 6 ) plane. As can be seen, a significant 
amount of parameter space is forbidden for the theory to be valid until the Planck scale. 

( 2 ) Similar to the previous section, the upper and lower bounds on A^ are placed from the requirements 
of perturbative unitarity and vacuum stability conditions respectively, as shown in Fig. 5. In this 
region too the lower limit of A^ has dependence on A 2 . The same happens for A 3 , as depicted in 
Fig. 6 . 

(3) It is to be noted that the theoretically allowed parameter space in this high DM mass region is 
fully consistent with the Higgs to diphoton LHC data. Since, a heavy H ± which naturally arises 
in this region, does not cause any significant deviation in diphoton signal strength for the SM-like 
scalar. This occurs even with a large positive A 3 (within the bounds shown in Fig. 6 ). In principle, 
a large negative A 3 could modify /i 77 unacceptably. However, as Fig. 6 shows, such values are 
inconsistent with the aforementioned theoretical constraints. 

(4) It is worth noting here that the parameter space valid until the Planck scale and corresponding to 
M — 110 TeV shrinks significantly for M = 10 9 TeV, as can be seen in Fig. 5 and Fig. 6 . This is 
because, the significant part of the parameter space is discarded by the stability condition at large 
M where the neutrino Yukawa couplings becomes large 0(1O -1 ). Such large Yukawa coupling 
contributes to the beta function of A 2 as shown in Eq. 24b through the terms +A 2 ^ and — y that 
either makes A 2 perturbative in some cases or A 2 negative and the vacuum unstable in other. We 
elaborately discuss this feature below with some benchmark points. 


We have selected two benchmark points(BPl and BP2) as samples out of the allowed regions consistent 
with the relic density constraints. These points demonstrate how the different theoretical constraints 
can truncate the scale of validity of this scenario for two different values of right-handed neutrino 
masses M — 110 TeV and M = 10 9 TeV respectively. The parameter values are listed in Table. 1. 
BP1 and BP2 yield Ddm ^ 2 = 0.1151 and 0.1207 respectively, which is within with the allowed range 
of relic density. For M — 110 TeV, BP1 ensures a stable vacuum till the Planck scale (Fig. 7a). It is 
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M h (GeV) M h (GeV) 


Figure 5: Region(s) allowed in the Mh~ Al plane obeying the various constraints for M — 110 TeV (left 
panel) and M — 10 9 TeV (right panel) with three different choices of A uv and two values of A 2 (upper 
and lower panel). The full region (marked by RC + SI) (magenta) is allowed by the DM constraints 
alone. The overlapped regions labeled by A (red), B (cyan) and C (green) are consistent with the 
theoretical constraints up to A uv — 10 6 , 10 16 and 10 19 GeV, respectively. 
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M h * (GeV) M h * (GeV) 


Figure 6 : Region(s) allowed in the M^- A 3 plane obeying the various constraints for M — 110 TeV (left 
panel) and M — 10 9 TeV (right panel) with three different choices of A uv and two values of A 2 (upper 
and lower panel). The full region (marked by RC + SI) (magenta) is allowed by the DM constraints 
alone. The overlapped regions labeled by A (red), B (cyan) and C (green) are consistent with the 
theoretical constraints up to A uv — 10 6 , 10 16 and 10 19 GeV, respectively. 
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BP 

M h 

M h ± 

M a 

Al 

A 2 

BP1 

850.0 GeV 

854.0 GeV 

858.0 GeV 

0.02 

0.1 

BP2 

710.0 GeV 

712.0 GeV 

711.0 GeV 

0.11 

0.1 


Table 1: Benchmark values (BP) of parameters affecting the RG evolution of the quartic couplings. 
For each BP, two values of M, namely, 110 TeV and 10 9 TeV, have been used. 


also consistent throughout with perturbativity and unitarity. However, one has y u = 0.168 at M — 10 9 
TeV. For this value, the term 0(\2y%) has a dominant role in the RG evolution and A 2 starts increasing 
rapidly from 10 9 TeV on wards (Fig. 7b). 




(a) (b) 

Figure 7: RG running of different scalar quartic couplings corresponding to BP1. The solid, dashed, 
dashed dotted and dotted lines denote the evolution curves of the stability conditions vscl, vsc2, vsc3 
and vsc4 respectively. 


For M — 110 TeV, BP2 exhibits similar features in the evolution trajectory as in BP1 (Fig. 8a). For 
M = 10 9 TeV, the Yukawa coupling y„ starts with an initial value around 0.51. This is accounted for by 
the very small mass splitting, of the order of a GeV, between H and A. Thus, the dominant contribution 
from the RH neutrinos comes from the 0 (y ^) term that causes A 2 to drop sharply (Fig. 8b). Hence, in 
BP2, vacuum stability is destroyed shortly after 10 9 TeV as the condition A 2 > 0 gets violated. This 
particular feature can only be witnessed in the case of closely spaced Mh and M^, which is the primary 
requirement to satisfy the relic density constraints discussed before. This completes the explanation 
of how the allowed area can shrink due to different theoretical constraints for M = 10 9 TeV. 
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(a) (b) 

Figure 8: Same as Fig. 7 but for the benchmark point BP2. 

7 Summary and Conclusions 


We have examined the high-scale validity of a scenario that (a) offers a scalar dark matter, (b) radia- 
tively generates Majorana masses for neutrinos, and (c) is responsible for leptogenesis. For this, we 
extend the SM fields with one additional inert Higgs doublet field ($ 2 ) and three right handed neutrinos 
(Ni). These new particles are odd under a discrete Z 2 symmetry, while all the SM particles are even. 
Because of this discrete symmetry, $2 does not acquire any vacuum expectation value (vev) and has no 
tree-level couplings to fermions. In this scenario, one has five physical scalars (h, H , A, where, h 
is denoted as the SM like Higgs boson with a mass of 125 GeV. The lightest state between H and A is 
the dark matter candidate due to built in Z 2 symmetry. In our analysis we have assumed H to be the 
dark matter candidate. Neutrino masses are generated at the one-loop level. The neutrino masses and 
mixing angles are determined in terms of Yukawa couplings (y^), new Higgs particle masses (M#, Ma) 
and three heavy Majorana masses (Mg 2 , 3 ). In our numerical analysis we have assumed M\ is mass 
of the lightest state and considered two values, namely, M\ = M = 110 TeV and 10 9 TeV. These two 
mass scales are consistent with Leptogenesis . For simplicity, in our analysis, we have considered only 
one diagonal Yukawa coupling and to determine the value of this coupling we have scanned over M# 
and Ma for a given value of M, by keeping M v ~ 0(0.1 eV). 

The parameter space of this model is determined in terms of additional Higgs boson masses, 
Mh,Ma,M h ±, one quartic coupling A 2 and a set of quartic coupling combinations, A l. While scan¬ 
ning the parameter space of this model, we have imposed the LEP bound on scalar masses, Mh,Ma 
and M h ± . 

As far as the dark matter constraints are concerned, we have used the relic density limits obtained at 
3 a uncertainty from the PLANCK experiment and the direct detection cross-section limit from the 
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LUX experiment. Finally, for Mh < Mh/ 2, which would lead to large invisible decay width of the SM 
like Higgs boson, we demanded that the corresponding branching ratio is less than 20% as obtained 
from the model independent Higgs precision analysis. 

With these set of constraints in hand, we have then scanned the IDM parameter space for two different 
ranges of dark matter masses, 50 GeV < Mh < 90 GeV and Mh > 500 GeV. It should be noted that 
with Mh fixed at 125 GeV, A 2 , Au, Mh , Ma and M H ± determined all the remaining quartic couplings. 
We have used these quartic couplings as the electroweak boundary condition by setting the starting RG 
running scale Q = M t and evolved these couplings up to the scale Auv, where this scenario remained 
consistent. In the RG evolution of these quartic couplings, the neutrino Yukawa couplings started 
playing its role from the right handed neutrino mass scale Q — M on wards. Following are the salient 
features of this model that our analysis brings out. 


• In our study we have explicitly demonstrated that at the low DM mass region, the vacuum 
stability, perturbativity and unitarity constraints put stringent limits on the low-energy value 
of the coupling Au and the Z 2 ~odd scalar masses M h ±,Ma- These bounds strongly depend 
upon the scale up to which the theory is valid and the right-handed neutrino mass scale. It 
is interesting to note that all the parameter space allowed by the DM relic density and direct 
detection constraints lies well within the region allowed by the theoretical constraints valid up 
to the Planck scale. However, once we have imposed the condition that the Higgs to diphoton 
signal strength (/i 77 ) should lie within 2 a of the weighted average value of the ATLAS and CMS 
data, the allowed parameter space further squeezed. 

• The scenario with high DM mass region (Mh > 500 GeV) is significantly different from that 
of the low DM mass region. In this case, the DM being very heavy, the constraints from direct 
detection is rather insignificant. On the other hand, the relic density constraint is ensured by a 
degenerate Z 2 ~odd scalars (AM ~ 10 GeV). As a result of these, a large part of the parameter 
space in Al — Mh and A 3 — M H ± planes remain unconstrained. 

• However, the study of the high scale validity of high DM mass region has interesting consequences. 
The DM-allowed region reduces substantially after imposing the theoretical constraints and this 
reduction is strongly dependent on the cut-off scale (Auv). The effect of neutrino Yukawa 
couplings in the RG evolution of the different quartic couplings are also prominent in this case. 
We have found that with the increase in the right handed neutrino mass scale M, the neutrino 
Yukawa coupling (y u ) also increases, which in turn further reduces the allowed parameter space 
by either destabilizing the vacuum or violating the perturbativity bound. There is nonetheless 
a clearly identifiable region of the parameter space, which keeps the model valid all the way up 
to the Planck scale. This scenario is consistent with the measured Higgs-to-diphoton rates as 
measured during the 8 TeV run of the LHC. 
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Appendix 


A One-loop Renormalization group (RG) equations 


The RG equations for the gauge couplings, for this model, are given by [42], 

2 dds 


167 T 


dt 


= -7gi, 


^ 2 ! - - 3 * 3 - 
i6 * 2 % = 79 ' 3 - 

Here g', g and g s denote the U(l), SU(2)^ and SU(3) C gauge couplings respectively. 
The quartic couplings A i (i — 1,..., 5) evolve according to, 

1671-2 ^ = 12A2 + 4A2 + 4As A4 + 2A2 + 2A 5 + \ ( 3 ^ 4 + 9 ,a + 2 g 2 g' 2 ) 

-Xi(9g 2 + 3 g' 2 - 12y t 2 - 12 yl - 4y 2 ) - 12 yf , 

167f2 ^ = 12A i + 4A i + 4A 3 A 4 + 2Xj + 2A| 

+ |(3fl» 4 + g /4 + 2 g 2 g' 2 ) - 3A 2 (3 g 2 + g' 2 - i y 2 v ) - 4yl , 

2 d\ 3 _ , o\.\ , A \2 , o \2 , o \ 2 , 3 /g^,4 , „/4 


.2 _/ 2 \ 


167 T — (Ai + A 2 ) ( 6 A 3 + 2 A 4 ) + 4A 3 + 2A 4 + 2A 5 + -(3g + g —2 g gj 


I 67 r 


; dA 4 

dt 


)d\$ 


- a 3 (9 g 2 + 2>g' 2 - 6 y t 2 - 6 yl - 2 y 2 -2yl), 

= 2 (Ai + A 2 ) A 4 + 8 A 3 A 4 + 4 A 4 + 8 A 5 + 3 g 2 g' 2 

-A 4 (9 g 2 + 3 g' 2 - Gy 2 - Gy 2 - 2 y 2 - 2 yl ), 


(23a) 

(23b) 

(23c) 


(24a) 

(24b) 

(24c) 

(24d) 


I 67 t 2 — = (2Ai + 2A 2 + 8 A 3 + 12A 4 ) A 5 — \§(9g 2 + 3 </ 2 — Gy£ — 2y 2 — Gy 2 — 2y 2 ), (24e) 


dt 

For the Yukawa couplings the corresponding set of RG equations are, 

m2 lt = » (-«»? - -/ - ^a' 2 + \vt + + \y> 

= » (-89? - -/ ~ i»' 2 + \vt + yi +1 yl 

= Vt (“|S 2 - + +3»? + 1j5 + ^9? 


(25a) 

(25b) 

(25c) 

(25d) 
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